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ABSTRACT 

A  simple  linear  programming  formulation  is  given  for  finding  an  A*/-ball 
with  largest  radius  contained  in  a  polyhedral  set  defined  by  m  linear 
inequalities.  The  linear  program  also  has  m  linear  constraints  similar  to 
those  defining  the  set.  It  is  shown  that  finding  the  largest  ball  is  not  much 
more  difficult  than  finding  a  feasible  point.  When  the  center  of  ball  is 
fixed,  the  largest  radius  is  easily  obtained  as  the  smallest  of  m  ratios. 

The  results  can  be  extended  to  balls  defined  by  other  norms  such  as  elliptic 
norms . 
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SIGNIFICANCE  AND  EXPLANATION 


Finding  the  largest  ball  or  hypercube  in  a  polyhedral  set  has  many 
applications  in  operations  research.  This  work  gives  a  simple  linear 
programming  formulation  for  solving  the  problem.  The  effort  needed  to  solve 
the  linear  program  is  almost  the  same  as  that  for  finding  a  feasible  point  in 
the  given  polyhedral  set.  Hence  the  result  of  this  work  can  be  considered 
optimal . 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


B(y,Yip)  is  an  ordinary  ball  for  p  =  2.  For  p  =  “,  it  is  a  hypercube.  Intuitively 
if  F  contains  no  interior  points  the  largest  ball  will  have  Y  “  0,  otherwise  Y  >  0 
(including  the  case  Y  *  +*•).  Since  the  iP-ball  is  convex,  B(y,Y;p)  £_  F  if  and  only 
if  F  contains  all  the  extreme  points  of  8<y,Y,-p).  In  the  case  that  p  »  +",  B(y, yip) 


extreme  pc  ints. 


y  +  Yzk,  k  =  1,2,.. .,2n 


where  zk  =  [±1 ,±  1 , . . . ,±1 J T>  Therefore  the  problem  of  finding  the  largest  ball  can  be 
formulated  as 

(2)  max  Y  subject  to  aT(y  +  Yzk)  i  bi#  i  =  1,...,m,  k  =  1,...,2n 

(y,f) 
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Although  this  is  a  linear  program,  it  contains  2  *rr  constraints.  Hence  (2)  is 
practically  intractable  even  when  n  is  as  small  as  20. 

In  this  paper,  we  give  a  linear  program  formulation  with  only  m  constraints.  The 
linear  program  (see  (5))  is  very  similar  to  and  no  more  difficult  than  the  following 

formulation  for  finding  a  feasible  point  in  F. 

T 

max  y  subject  to  aty  — y  £  b^  i  »  1,...,m 

(y.Y ) 

Hence  the  problem  can  be  solved  by  efficient  algorithms  such  as  Dantzig's  simplex  method 
or,  if  the  problem  is  large  and  sparse,  Mangasarian' s  SOR  method  [6).  This  result  also 
shows  that  theoretically  the  problem  is  polynomial-time  solvable  (3),  (4).  It  is 
interesting  to  note  that  finding  the  smallest  ball  containing  F  is  much  more  difficult. 
Depending  on  the  norm  used,  it  can  be  NP-complete  (7) ,  which  means  that  if  one  can  solve  it 
in  polynomial  time  then  he  can  solve  also  in  polynomial  time  hundreds  of  those  intractable 
problems  such  as  traveling  salesman  problem  (2)  or  non-convex  linear  complementarity 
problems  (1).  These  problems  are  considered  intractable  because  as  it  is  widely  believed, 
but  not  proven,  that  no  polynomial-time  algorithms  exist  for  solving  them. 

2.  LP  Formulation 

The  problem  can  be  written  as 

(3)  maximize  Y  subject  to  B(y,y,p)  F 

(y.T ) 

For  i  “  1,2,...,m,  define  the  function  g^(y»Y»p)  by 

(4)  g1(y.Y»p)  s“  max  a^x  -  subject  to  x  e  B(y,y»p) 

x 

It  is  easy  to  see  that  the  constraint  in  (3)  can  be  replaced  by  m  constraints 

g1(y»t»p)  £  0  i  =  1,2, ...,m. 

Lemma  1  :  B(y,y;p)  c  F  if  and  only  if  9^(y»Y;p)  £  0  for  i  “  1,2, ...,m. 

Proof :  B(y,Y;p)  c  f 


aT(  x-y ) 


a 

1  i3 


q-1 


<Z  la 


1 3 


q,  1-1/p 


=  Y 


( Z  |  a  .  ,|q 

ij' 


) 1/q  =  Y 


|a.  | 

1 3  q 


It  follows  that  (3)  is  equivalent  to 

(5)  maximize  Y  subject  to  a?y  +  Y*|a.  |  ^  b^,  i  =  1,2,  . .  .  ,m 

(y »Y )  1  q 

We  summarize  the  results  in  the  followsing: 


□ 


Theorem  1 .  For  any  given  a^,  b^ ,  i  =  1,2,...,m 

(i)  The  linear  program  (5)  is  feasible. 

*  * 

(ii)  Assume  that  the  linear  program  (5)  is  bounded,  and  that  (y  )  is  an 

* 

optimal  solution  with  Y  <  +°*.  Then 

(a)  Y*  <  0  if  and  only  if  F  » 

(b)  Y*  =  0  if  and  only  if  F  ^  $  but  F  has  empty  interior. 

(c)  Y*  >  0  if  and  only  if  F°  /  ♦  and  B(y*,Y*;p)  is  an  fP-ball  contained 

in  F  with  the  largest  radius. 

Proof .  Since  (5)  is  equivalent  to  (3)  by  Lemma  1,  2.  The  theorem  follows  by  the  following 
observations. 

(i)  (0,-K)  is  feasible  for  K  sufficiently  large. 

(ii)  (y,0 )  is  feasible  for  all  y  e  F. 

(iii)  F  has  non-empty  interior  if  and  only  if  F  contains  a  ball  B(y,Y;p)  with 
Y  >  0. 

□ 

* 

Remarks  (i)  If  Y  =  +*  then  F  has  unbounded  interior.  But  the  reverse  is  not  true. 
For  example,  let 

F  =  {(x1,x2)|x1  -  x2  <  1,  -x1  +  x2  £  0} 


The  same  example  also  shows  that  y  is  not  unique. 


(ii)  Given  y  e  F,  the  largest  ball  contained  in  F  and  centered  at  y  can  be 


(6) 


found  by  solving 

max  Y  s.t.  Y* I I _  <  bj.  ~  a^y,  *■ 

V  X  J  " 


which  can  be  solved  explicitly,  namely. 


* 

Y 


min 

Ki<m 


.  T 
b .  -a  y 
l  l 


1, . .  .  ,m 


3.  General  Norms 

From  the  derivation  of  (5)  in  the  previous  section,  it  is  clear  that  the  results  can 
be  generalized  to  other  norms.  That  is,  find  a  largest  ball  B(y,y)  =  {x  |  |x-y|  i  y} 
in  F  is  equivalent  to 


(7) 


max  y  s.t.  a*y  +  y*|a.|  &  b^,  i  =  1,...,m 

(y.Y)  1 


where  | |  is  the  dual  norm  of  | “ | ,  namely 


(8) 


l^Z  s.t.  I  z  1  <  1 


Of  course,  to  make  (7)  useful  computationally,  we  need  to  be  able  to  solve  (8),  as  in 
Section  2.  For  example,  if  |z|  :=  (zTA where  A  is  a  symmetric  positive  definite 

matrix,  then  by  Kuhn-Tucker  Theorem  (5]  we  have 

i  i  *  T  T 

a.  =  max  a.z  s.t.  z  Az  <  1 

z  1 


(9) 


,  T-1  ,1/2 

=  ( a  .  A  a  .  ) 
i  i 


Note  that  in  this  case,  a  "ball”  is  in  fact  an  ellipsoid.  Hence  one  can  find  the  largest 
ellipsoid  (with  a  given  shape  defined  by  A)  in  F  by  solving  a  linear  program.  When  the 
center  is  fixed,  the  largest  ellipsoid  can  easily  be  found  by  (6)  in  which  lajjq  is 
replaced  by  (a^A 
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